To APPEAR IN THE ASTROPHYSICAL JOURNAL LETTERS 
Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



O 

o 

(N 

I 

o 



> 
o 

o 
a^ 
o 



NO NEUTRON STAR COMPANION TO THE LOWEST MASS SDSS WHITE DWARF 

Marcel A. Agueros^'*, Craig Heinke^, Fernando Camilo\ Mukremin Kilic'' '\ Scott F. Anderson'', Paulo 
Freire"', Scot J. Kleinman"*, James W. Liebert^, Nicole M. Silvestri'' 

To appear in the Astrophysical Journal Letters 

ABSTRACT 

SDSS J091709. 55+463821. 8 (hereafter J0917+4638) is the lowest surface gravity white dwarf (WD) 
currently known, with log g = 5.55 ± 0.05 (M « 0.17M,t.: iKilic et al.|[200 7a,bV Such low -mass white 
dwarfs (LMWDs) are believed to originate in binaries that evolve into WD/WD or WD/neutron star 
(NS) systems. An optical search for J0917+ 4638's companion showed that it must be a compact object 
with a mass > 0.28 Mq (jKihc et al.ll2007bD . Here we report on Green Bank Telescope 820 MHz and 
XMM-Newton X-ray observations of J0917-I-4638 intended to uncover a potential NS companion to 
the LMWD. No convincing pulsar signal is detected in our radio data. Our X-ray observation also 
failed to detect X-ray emission from JG91 7-1-4638 's companion, while we would have detected any of 
the millisecond radio pulsars in 47 Tuc. We conclude that the companion is almost certainly another 
WD. 

Subject headings: white dwarfs — stars: individual (SDSS J091709. 55-1-463821. 8) — pulsars: general 



1. INTRODUCTION 

Low- mass white dwarfs (LMWDs), generally defined 
as having M < 0.45 Mq, make up a small but highly in- 
teresting subset of white dwarfs (\yDs). Using the Palo- 
mar Green Survey, iLiebert et al.l ()2005f ) estimated that 
the formation rate of LMWDs is 0.4 x 10"^^ pc^^ yr"^, 
meaning that they make up only '^10% of the population 
of the most commonly observed WDs, hydrogen atmo- 
sphere DAs. But it is their presumed evolutionary histo- 
ries that make LMWDs truly intriguing. The youngest 
WDs in the oldest globul ar clusters in the M ilky Way 
have masses of ~ 0.5 Mq (|Hansen et al.l 120071 ). implying 
that lower mass WDs undergo significant mass loss as 
they form. The preferred scenario is that these WDs 
form in a tight binary whose evolution includes a phase 
of mass transfer. As a result, much of the WD progen- 
itor's envelope is removed, preventing a helium flash in 
it s core, and producing a low-mass, helium-corc WD. 

iBrown et al.l (|200l identified J09 17+4638 as a DA 
WD in their hyper-velocity star survey of photometri- 
cally selected B-st ar candidates. Deta iled model atmo- 
sphere analyses by IKilic et ahl (|2007al lbl) showed that it 
has Teff = 11,855 K, log g = 5.55, and M w O.17M0. 
The lack of evidence of a companion in the optical pho- 
tometry forces any main-sequence companion to have 
M < 0.1 Mq, ruling out a low- mass main-sequence star 
companion. Radial velocity monitoring uncovered varia- 
tions with a period of 7.6 hr, implying t hat the mass of 
the companion is > 0.28 Mq (jKilic et al1l2007bf) . 
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What is the nature of this compani on? While LMWD s 
are found in WD/WD systems fe.g.. lMarsh et al.|[l995D . 
most known LMWDs are found as companions to neu- 
tron stars (NSs), and specifi cally to NSs "recy cled" as 
millisecond pulsars (MSPs; iPanei et all l2007t i. Most 
field radio pulsars in binary systems are MSPs, where 
a middle-aged, radio-quiet NS has been reactivated as 
a pulsar via accretion from its companion. The MSP 
companions are generally thought to be LMWDs with 
M = 0.1 — 0.4 Mq, although they are often too faint for 
opti cal spectroscopy to confir m that they are LMWDs 
(see Ivan Kerkwiik et al.ll2005f) . Still, a third of the ~50 
MSP companions discovered outside of globular clusters 
have M <, 0.2 Mq, assuming the systems have a median 
inclination of 60° (Manchester et al. 2005b). 

While simulations designed to identify the evolution- 
ary pathways that produce LMWD /MSP systems do not 
generally predict many systerns wit h Porb much shorter 
than a day (e.g.. iNelson et al]|2004f ). and while the sys- 
tem's mass function implies that the probabilit y that 
J0917+ 4638 has a WD companion is 89% ( Kilic et all 
I2007bf) . a NS (or black hole) companion to this LMWD 
cannot be ruled out with the current optical observa- 
tions. In addition, for the currently known sample of 
WD/WD systems where both WD masses have been 
measured, the mass ratio is typically about unity (see 
iNelemans et al.l (20051 and references therein), while the 
ratio for the J0917+4638 binary system is < 0.61. 

Because of the connections between LMWDs and 
MSPs, we used the Green Bank Telescope (GBT) to 
search for a putative pulsar companion to J09 17+4638, 
and report here on these observations (Section [2]). We 
also report on an XMM-Newton X-ray Observatory ob- 
servation of this LMWD (Section[3l) . Blackbody emission 
from a putative NS companion to the LMWD should 
be gravitationally bent, allowing us to detect the NS in 
X rays even if it were radio-q uiet or if its pulsar beam 
were missmg our line of sight (jBeloborodovl l2002f ) . We 
are specifically motivated by the X-ray de tection of all 
known MSPs in the globular cluster 47 Tuc (iHeinke et al.l 
[2005HBogd anov et al]l2006f ). allowing predictions of the 
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TABLE 1 

SDSS J0917+4638: Properties and Observations 



SDSS g Teff MwD -Porb Dist. b DM GBT XMM 

(mag) (K) {Mq) (hr) (kpc) (°) (cm-^ pc) Int. (s) Int. (s) 

18.77 ill 0.02 11855 0.17 7.594 ± 0.002 2^3 +44.0 80 13300 23418 

Note. — The g (PSF) magnitude is from SDSS Data Release 7 l|Ab azajian et al.ir2009l ). The distance and orbital period are from [Kilic et all 
<|200 7b.^ ■ The listed DM is the maximum value used when searching f or pulsations; it corresponds approximately to twice the maximum value 
obtained in tlie direction of J0917+4638 with the lCordes &: Laziol I I2002I) model. 



X-ray emission of other MSPs. We choose this sample 
of MSPs for comparison in part because the distance to 
globular clu sters such as 47 Tuc (4.5 kpc, 2003 update of 
IHarrislHOOl are better known than the distances to most 
MSPs. We discuss the significance of our nondetections 
and conclude in Section [H 

2. GREEN BANK TELESCOPE OBSERVATIONS 

J0917+4638 was observed with the GBT on 2007 
November 30. The observing set-up and data re- 
ductio n were the same as described in lAgiieros et alj 
()2009t ) At 820 MHz, th e Berkeley-Caltech Pulsar Ma- 
chine (|Backer et al.lll997f ) provided 48 MHz of bandwidth 
split into 96 spectral channels; total power samples for 
each channel were recorded every 72 fxs. The total ob- 
serving time was 13,300 s (3.7 hr). We used standard 
pulsar search techniques as implemented in the PRESTO 
software package (Ransom 2001). We calculated the 
maximum dispersion measure ( DM) expected in t he di - 
rection of J09 17-1-4638 using the ICordes fc Laziol (j2002l ) 
model for the distribution of free electrons in the Galaxy. 
To account for uncertainties, we dedispersed the data up 
to a DM limit twice that obtained from the model, i.e., 
DM = 80 cm-3 pc. 

No convincing pulsar signal was detected in our data. 
Below we discuss the limitations of our search. 

2.1. Acceleration Sensitivity 

The orbital motion of a putative pulsar companion to 
J0917-I-4638 could significantly affect its ap parent spin 
period. Based on radial velocity monitoring. iKilic et al.l 
(|2007bf) found that J09 17-1-4638 is in an orbit with a pe- 
riod 7.6 hr. Assuming that the LMWD companion is 
a 1.4 Mq NS, this implies that the maximum orbital 
acceleration is on the order of 100 m s~^, which is sig- 
nificantly larger than what is typically seen in these sys- 
tems (for 90% of known puls ars the maximum orbital 
acceleration is < |25| m s~^; [Manchester et all l2005bl : 
Ivan Leeuwen et ani2007t l. 

Our integration time represents nearly half of the bi- 
nary orbital period. As a result, the assumption of a con- 
stant apparent acceleration built into PRESTO breaks 
down. We therefore divided our GBT data into 14 sep- 
arate 900 s integrations (each representing ^ 3% of an 
orbit) and one 700 s integration and conducted searches 
for pulsations separately in each of these partial obser- 
vations."'^'^ This extended our search sensitivity to accel- 
erations on the order of several hundred m s~^, but as 
detailed in the following section, reduced our luminosity 
sensitivity. None of these searches uncovered a convinc- 
ing pulsar signal. 

We also conducted a search of the entire 3.7 hr integration, 
which unsurprisingly did not return any good candidate pulsar 
signals. 



2.2. Luminosity Sensitivity 

We use the standard modifications to the radiometer 
equation to calculate the minimum detectable period- 
averaged flux density for our searches. We consider a 
pulsar duty cycle of 20% (typical of MSPs). At 820 MHz, 
the GBT gain is 2K Jy^^ and the system temperature is 
25 K. The sky temperature at this frequency and a Galac- 
tic latitude of 6 = -1-44° only adds a few K to the overall 
temperature. We consider an effective threshold signal- 
to-noise ratio of 10. For tint — 900 s, the sensitivity limit 
for a long period pulsar at the beam center is ~0.26 mJy. 

Pulsar luminosities are often measured at 1400 MHz; 
using a typical spectral index of —1.7, the limiting sensi- 
tivity at that frequency is S'i4oo ~ 0.10 mJy when search- 
ing the 900 s integrations. For an MSP period of 3 ms, 
our sensitivity at 1400 MHz was roughly 0.14 mJy for 
each integration, and it quickly degraded for shorter pe- 
riods; it was 10 X worse for 1 ms. 

The distance to J 09 17-1-4638 is estimated to be 2.3 kpc 
(|Kilic et al.l l2007bl ) . implying that our L1400 = 5*1400^^ 
hmits for 3 ms periods are « 0.7 mJy kpc^ for each 900 s 
integration. According to the ATNF's pulsar catalog^^ 
(|Manchester et al.|[2005ah . of the 50 MSPs (periods < 
25 ms) outside of globular clusters and with measured 
luminosities, 64% have L1400 > 0.7 mJy kpc^. We would 
therefore expect our search to detect roughly two-thirds 
of the known MSPs were one orbiting J0917-I-4638 and 
beaming radio waves toward the Earth. 

We note that J0917H-4638 falls within the FIRST foot- 
print and is not detected in that 1.4 GH z survey, for 
which the sensitivity limit is roughly 1 mJy (|Becker et al.l 
fT995h . 

3. XMM-NEWTON OBSERVATION 
3.1. Motivation 

MSP radio beaming fractions are < 100%, and as a 
result, some MSPs have not yet been detected in the 
radio in binary systems where there is strong evidence 
for their presence (e.g. , the companio n to the young pul- 
sar PSR J1906-H0746: iLorimer et 31112006^ . Given that 
the NS blackbody emission is gravitationally bent, al- 
lowing us to view > 75% of the NS surfaces in X rays 
(|Beloborodovl l2002^. sufficiently deep X-ray observations 
a re virtually guarante ed to detect these MSPs. 

iHeinke et al.l (|2005( l found no correlation between the 
X-ray and radio luminosities of MSPs in 47 Tuc, as ex- 
pected due to the differing nature and spatial location 
of the X-ray and radio emission, and found that all 
MSPs in 47 Tuc^^ have X-ray luminosities ranging be- 
tween Ljf (0.5 - 6 keV) = 2 x 10^" and 2 x lO^^ erg s'^ 

http: //www. atnf . csiro . au/research/pulsar/psrcat/. 

Only 15 MSPs in 47 Tuc have published locations farther than 
1" from other MSPs; two pairs of MSPs that are closer cannot 



NO NS COMPANION TO J0917+4638 



3 



iBogdanov et al.l (|2006f ) showed that the X-ray spectra of 
the MSPs in 47 Tuc are typically dominated by thermal 
blackbody-like emission from the NS surface around the 
polar caps, with temperature 1 — 3 x 10^ K. This X-ray 
emission is sometimes overwhelmed by additional non- 
thermal X rays that are e ither magnetospheric o r due 
to an intra-binary shock. IBogdanov et al.l (|2006D also 
showed that there are no clear systematic differences be- 
tween the X-ray properties of MSPs in 47 Tuc and in the 
Galactic field. Thus, we requested an XMM observation 
capable of detecting any of the known MSPs in 47 Tuc, 
were they located at the distance of J09 17-1-4638. 

J09 17-1-4638 had not previously been obs erved in the 
X-ray since the ROSAT All-Sky Survey (jVoges et all 
[199112000), where it was not detected (unsurprisingly, 
considering the short exposure time). 

3.2. X-Ray Data Analysis 

We observed J0917-h4638 on 2008 May 7 for 17 ks (Ob- 
sID 0553440101) with XMM s EPIC camera, con sisting 
of two MOS CCD detectors (iTurner et al.| [200lD and a 
pn CCD detector (jStriider et a l.l l2001f ). All data were 
reduced using FTOOLS" and SAS version 8.0.0." We 
excluded times of soft proton background flaring, when 
the pn camera's count rate exceeded 25 (0.2 — 10 keV) 
counts s~^, or when the MOS cameras exceeded 7 or 8 
(0.2-10 keV) counts s~i for the MOSl or M0S2 cameras 
respectively. This left 8.9 ks of good data from the pn de- 
tectors, and 11.3 ks from the MOS detectors. We filtered 
the events on pixel patterns (trying PATTERN<=1 or 
<=4 for pn and PATTERN< = 12 for MOS data), and 
for FLAG==0. We choose an energy range of 0.2 — 1.5 
keV to obtain optimal sensitivity to the soft blackbody 
emission expected from MSPs. 

No source is detected at or within 1' of the loca- 
tion of J0917-f4638 either with detection algorithms or 
by eye. We utilize our knowledge of the XMM point 
spread function^^ and absolute pointing accuracy (< 1"; 
iKirsch et al.ll200^ to determine an upper limit. For the 
pn, 50% of 1.5 keV photons are found within 8", and 80% 
within 20". For the MOS cameras, 50% of 1.5 keV pho- 
tons are recorded within 8", and 75% within 20". We 
find 3 counts within an 8" circle, or 20 counts within 
a 20" circle, in the combined image. This is consistent 
with a nondetection, as the expected background counts 
in these circles are 3.2 ± 0.2 and 19.8 ± 1.3 counts, re- 
spectively, as derived from nearby background regions. 

3.3. Comparison to MSPs in 47 Tuc 

We use the X-ray faintest MSP in 47 Tuc, 47 Tuc 
T, to calibrate our expectations for the detection of 
an MSP in J0917-h4638. 47 Tuc T has Lx(0.2 - 1.5 
keV) = 1.5 X 10^° erg s~^ an d a 134 eV blackbody spec- 
trum (jBogdanov et al.ll2006[ ). We use PIMMS^*' to deter- 
mine the expected EPIC count rates from 47 Tuc T were 
it located at 2.3 kpc (the distance to J0917-I-4638) behind 
an estimated Nh = 1.5 x 10^" cm~^ (jDickev &: LockmanI 



Il990f) . We expect 10.9 counts within 8", or 17.0 within 
20", accounting for the relevant encircled energy frac- 
tions, from such an MSP. Comparing the predicted 
counts with th e Poisson errors on the detected counts 
(|Gehrelsl 11986'. equation 7), we find that the number 
of counts within 20" is 3. Oct below expectations for the 
faintest known MSP in 47 Tuc, while the counts within 
8" are 3.5cr below those expectations. 47 Tuc T is the 
X-ray faintest of the 15 independently measured MSPs 
in 47 Tuc; the median X-ray luminosity is 2.1 x greater 
(jBogdanov et al.ll2006D . which is ruled out at 6.3cr con- 
fidence. Our nondetection is therefore strong evidence 
against the existence of an MSP in J09 17-1-4638. 

4. CONCLUSION 

We have searched for evidence of an MSP companion to 
the LMWD J0917-I-4638 through radio and X-ray obser- 
vations. Our radio search reaches a sensitivity sufficient 
to detect roughly two-thirds of the known MSPs, while 
our X-ray search is sensitive enough to detect any of the 
15 independently identified MSPs in 47 Tuc. Together, 
our nondetections provide strong evidence against the 
presence of an MSP in this system. Furthermore, since 
any NS companion to J09 17-1-4638 would presumably 
have been recycled through accretion from the LMWD, 
we rule out the presence of a NS in this system. Although 
a black hole companion is still conceivable (as the 7.6 hr 
orbital period would not induce current accretion and X- 
ray activity), such a companion is far less probable than 
a WD companion g iven both the system's mass function 
(jKilic et al.ll2007b[ l an d the stellar initial mass function 
for M > 1 Mq (e.e... IScalolll998f) . We conclude that 
J0917-h4638's more massive companion (M > 0.28 Mq) 
is almost certainly another WD. 

Roughly two dozen WD/ WD binaries are known and 
in te n such systems both W D masses have been measured 
(see iNelemans et al.ll2005L and references therein). The 
individual masses of WDs in these systems range between 
0.29 and 0.71 Mq; the median mass for those with mea- 
sured masses (and not just lower limits) is 0.43 Mq. The 
majority of these double WD systems have mass ratios 
near unity, which is contrary to wha t is expected from 
standard population synthesis models (jNelemans fc Touti 
l2005f ). This has been used to argue that energy bal- 
ance (a-formalism) , the standard prescription for com- 
mon envelope evolution, should be replaced by angu- 
lar momentum balanc e (^-algorithm: iMaxted et al]l2002t 
INelemans eraIll2005D. 

In particular. INelemans et al.l (|2005f ) found that the a- 
formalism cannot be used to describe the first phase of 
mass transfer for nine of the ten double WD systems in 
which both WD masses have been measured. The excep- 
tion is WD1704-I-481, which has a mass ratio = 0.7, simi- 
lar to the expected ratio from the a-formalism. The mass 
ratio for the J0917-)-4638 binary system is < 0.61. Recent 
observations of another LMWD, LP400— 22, showed 



be conclusively resolved l|Bogdanov et aLll2006l '). However, the flux 
from each pair is consistent with expectations from the other MSPs. 

http: //heasarc .gsfc.nasa. gov/docs/software/f tools/ ftools_menu.html. 
http://xiimi.vilspa.esa.es. 
15 XMM-Newton User's Handbook, 

http : //xmm. esac . esa. int/external/ximn_user_support /document at ion/uhb/ index .html. 
1® http://asc.harvard.edu/toolkit/pimms . jsp. 
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that it is in a binary with a mass ratio < 0.46 (|Kilic et al.l 
l2009f ). The mass ratios of these systems imply that the 
Qf-formahsm may explain at least some of the WD/ WD 
binaries. Determining the mass ratios of the other SDSS 
LMWD systems for which the nature of the companion 
is currently unknown will be important in understanding 
the role of energy versus momentum balance in recon- 
structing common envelop evolution. 
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